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Abstract 

We report the numerical calculations of the static vortex structure and criti- 
cal curves in exponentially shaped long Josephson junctions for in-line and overlap 
geometries. Each solution of the corresponding boundary value problem is associ- 
ated with the Sturm-Liouville problem whose minimal eigenvalue allows to make a 
conclusion about the stability of the vortex. The change in width of the junction 
leads to the renormalization of the magnetic flux in comparison to the case of a 
linear one-dimensional model. We study the influence of the model's parameters 
and, particularly, the shape parameter on the stability of the states of the magnetic 
flux. We compare the vortex structure and critical curves for the in-line and overlap 
geometries. Our numerically constructed critical curve of the Josephson junction 
matches well with the experimental one. 

Introduction 

The Physical properties of the magnetic flux in Josephson junctions are the base of modern 
superconducting electronics. The shape of the junction has an essential influence on its 
properties [T]. One of the very interesting objects is the exponentially shaped Josephson 
junction (ESJJ). The unidirectional flow of the magnetic flux quanta in the long Josephson 
junction provides a means for tunable oscillators at the frequencies above 100 GHz with 
power output of a few microwatts. The ESJJ decrease the spectral linewidths of the 
oscillator, which is large in junctions with rectangular shape. The exponential shape 
provides better impedance matching to an output load and allows to avoid a chaotic 
regime that exists for rectangular junctions [7. 

Josephson junctions with an exponentially varying width in the xy plane have been 
studied theoretically and experimentally in where the influence of the junction's 

shape on the current-voltage characteristics of the junctions was studied in detail. How- 
ever, the problems involving the determination of the stability regions of the static dis- 
tributions and the structure of the critical curves have not been adequately studied. The 
present paper is devoted to a study of those questions, which are important from the 
applied and theoretical points of view. 



1 Method 



The basic equation for the phase tf{t,x) in the ESJJ can be written in the form j2] ^ |1] 

if + ai^ — ip" + sirup + g(t,x) + 'J = . (1.1) 

Here x G (0, 1) and t > 0. The spatial coordinate x is normahzed to the Josephson 
penetration depth, / is the dimensionless junction's length, 7 is external current and 
time t is normalized to the plasma frequency . An overdot denotes differentiation with 
respect to time t, and a prime denotes differentiation with respect to the spatial coordinate 
X. The parameter a describes the dissipative effects. The term 

g{t,x) = a[y^'{t,x)-hB] (1.2) 

in the equation (|1.1|) describes the additional current in the exponentially shaped Joseph- 
son junctions with the shape parameter a > 0. As it is known P, by appropriate normal- 
ization we can interpret (p{t, x) as dimensionless magnetic flux along the junction. Then 
hs represents the dimensionless intensity of the external (boundary) magnetic field. 

In this paper we consider the static distributions fix) of the magnetic flux along the 
junction. So, the corresponding boundary value problems (EVP's) are posed as follows 



^/ for the in-line geometry 



—(p" + siinp + g{x) = 0, (1.3a) 
v'{0)-hB + h = 0, (1.3b) 
^'{l)-hB = 0, (1.3c) 



a/ for the overlap geometry 



-if" + simp + g{x) + = , (1.4a) 
^'{0)-hB = 0, (1.4b) 
ip'{l)-hB = 0, (1.4c) 

The solutions of the EVP's (|1.3p and (jl.4|) depend on the set of four parameters 
p = {/, a, Hb, 7} of the model, i.e. ip = ip{x,p). 

In order to study the stability of the solutions ip{x), following j5j we pose [H] the 
Sturm-Liouville problem (SLP) 

+ a'ip' + q{x)^ = Xip, xeiOJ), (1.5a) 
^'(0)=0, ^'(0 = 0, (1.5b) 

def 

Here A is the spectral parameter. The potential q{x) = cos(y9(x) corresponds to the 
concrete static solution (p{x). 

We note that the function q{x,p) is limited |g'(a;,p)| < 1 on the finite interval [0,/] 
for every p E V C M."^. Then the SLP (jl.5|) have [7j bounded below counted sequence 
of real eigenvalues (EV) {A„}, n = 0, 1, . . . Every EV A„ corresponds to an unique real 
eigenfunction (EF) ipnix), which satisfies the norm condition 

ipni^) dx = 1 . (1.6) 
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The number of zeroes of EF ipn{x) in the interval (0, /) is equal to n. Particularly, the EF 
ipo{x), which corresponds to the minimal EV Xmin = Aq does not have zeroes for x G (0, /). 

Because the solution ip{x, p) depends on the parameters p, the potential of SLP induced 
by this solution and the corresponding EV and EF depend on the parameters as well, i.e.. 
An = K{p) and ipn = i^n{x,p). 

The static solution is exponentially stable [5, 6^ in some subset of the parameter region 
V if the minimal EV satisfies Xminip) > 0. When Xminip) < the solution ip{x) is unstable. 
The points p eV lying in the hypersurface 



are the bifurcation points for the solution ip{x, p) . The values of the parameters that satisfy 
(ll.7|l are called critical (bifurcation) values for the solution ^p{x). The cross section of the 
surface ()1.7|1 by hyperplane corresponding to fixed values of two parameters determines 
the critical curve for another two parameters. From an experimental point of view, the 
most important are the critical curves "current-magnetic field" 



for fixed length / and shape parameter a. 

We have to note that mathematically it is possible to pose two different problems for 
the equations ()1.3|) - p.6p . 

Let '^{x) be a solution of BVP (jl.3|) or (jl.4j) corresponding to a fixed set of param- 
eters p. We can check its stability by means of SLP (|1.5|) . Let us suppose that 
Xminil, hs) > 0. In order to obtain a point from the critical curve (jl.Sp . we have to 
change for given 7 (or change 7 for given hs) until the condition p.Sj) is fulfilled. 

In order to calculate bifurcation points directly, we have to solve the system p.3p or 
()1.4|1 . and for given Xmin as a non-linear eigenvalue problem with respect to 

the couple {(p{x),ip{x)} and one of the parameters p |H] ^ [IDI- Such an approach has 
been successfully applied to different nonlinear physical problems [TT] . 

In this paper, we find the regions of stability for distributions of magnetic flux in ESJJ 
by numerical calculations. Let us note that the replacement of the problem for stability of 
the solutions of nonlinear operator equations by the eigenvalue problem for linear operator 
has the rigorous mathematical explanation ^21, [T^ - 

2 Results and discussion 

The curves of Amin(^B) for the Meisner distribution and the first few stable vortices in 
a Josephson junction in the in-line geometry for current 7 = and two values of shape 
parameter a = and a = 0.07 are demonstrated in Fig. ^ 

Each curve has two zeroes, the distance between which determines the stability region 
of the vortex upon variation of the magnetic field hs- The zeroes themselves are critical 
values of the field hs for zero current 7. It is well seen that the role of the shape parameter 
a is most important for small values of hs- In particular, when cr = the $^ vortex exists 
already for hs ~ 0.054. For a = 0.07, however, the existence region in respect to field Hb 
is considerably compressed, and the vortex exists starting for hs > 0.75. The compression 
of the Amin(^B) curves for different vortices decreases rapidly with increasing hs- 



Xminip) = 



(1.7) 



Xminil, hs) = , 



(1.8) 



3 




12 3 4 

Boundary Magnetic Field 



1.0 



0.9 



C 0.8 
BO 



W 



0.0 







Case: l = l,G 


= 0.07,7 = 


- 












M 
















T?f 

1 1 1 1 1 1 





1 2 3 4 5 

Boundary Magnetic Field 



Figure 1: The Aniin(^s) curves for different a 

The dependencies of the minimal eigenvalue of the SLP for the fluxon $^ on the shape 
parameter a in case of in-line geometry for hs = 1 and hs = 2 are shown in Fig. El It is 
seen that for fixed hs and 7 there exists a certain maximum value of a above which the 
distribution of loses its stability, i.e., a bifurcation of the vortex occurs upon a change 
in a. Large values of the magnetic field Hb correspond to large critical values of a. 
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Figure 2: Bifurcation upon a change in a for /i^ = 1 (left) and hs = 2 (right) in-line 
geometry 

The influence of the external magnetic field hs on the magnetic flux distribution ip '{x) 
in the in-line geometry for the fiuxon $^ for a = 0.07 was studied in 6J. For a certain 
value Hb = ~ 1.515, the maximum of the derivative ^'{x) is localized in the middle of 
the junction. For hs < hm the fiuxon is "expelled" to the end x = / by the "geometric" 
current g{x). If the length of the Josephson junction is sufficiently large, then a change 
of the current 7, equivalent to a change in magnetic field for the left end x = 0, turns 
out to have only a slight infiuence on the local maximum magnetic field in the junction. 
For (7 = the variation of the current causes the maximum of the magnetic field to shift 
to the right or left of the center of the Josephson junction, depending on the direction of 
the current. For values hs > hm the scheme Xm of the maximum of the field (p'{x) of the 
fiuxon $^ is shifted to the left from the center x = 1/2 toward the end x = 0. 

As we can see in Fig. |2l the stability region in case of in-line geometry is increased for 
hs = 1 {hB < hm) and decreased ioi Kb = 2 {Iib > hm) with the increase of the external 
current. 
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For the overlap geometry (see Fig. El) we have the opposite situation. The stabihty 
region is decreased for /i^ = 1 {hs < hm) and increased for hs = 2 {hs > hm) with the 
increase of the external current. 




Shape Parameter Shape Parameter 



Figure 3: Bifurcation upon a change in a for hs = 2 (left) and hs = 2 (right) in overlap 
geometry 

When 7 = the BVPs for the the in-line and overlap geometries coincide. In case of 
in-line geometry the increase of 7 leads to the decrease of /i^ at the end x = 0. The sign 
of g{x) depends on the value of hs and for hs > hm its sign is negative. When hs < hm, 
there is a region on x where g{x) > 0. In the first case, the EVP leads to the increase 
of the stability region with the increase of 7, in the second one - to its decrease. In the 
overlap geometry the boundary conditions do not change, but the additional term 7 in 
the EVP fll.4|) effectively changes the absolute value of g{x). For g{x) < {hs > hm) the 
stability region is increased with 7, because 7 decreases the absolute value of g{x). When 
hs < hm external current increases the interval in which g{x) > and this leads to the 
decrease of the stability regions as shown in Fig. IHl (left). 
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Figure 4: The Amin(7) curves for the in-line (left) and overlap (right) geometries 

Let us now consider the question of constructing by numerical means the critical 
current versus magnetic field relation, which is determined implicitly by Eq. (jl.8|) for 
each magnetic flux distribution in the Josephson junction. The importance of this problem 
stems from the possibility of measuring this relation experimentally El El El E] • 
We note that for different configurations of the magnetic flux in the Josephson junction 
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the values of the critical parameters (in particular, the critical current and magnetic field) 
can be substantially different. Therefore we have to distinguish the critical parameters 
for the specific distributions and for the Josephson junction as a whole. 

Figure El shows the Amin(7) curves for stable solutions of the BVP (jl.3|) in a field 
Hb = 1.7 for both in-line and overlap geometries. The distances between zeroes of the 
functions are the stability intervals of the corresponding distributions with respect to the 
current 7. The right and left zeroes of the function Ainin(7) are the positive and negative 
critical currents, respectively, of the distribution in the given field hs- Because of the 
asymmetry of the boundary conditions in in-line geometry for 7 > the critical current 
of the Meissner distribution [which we denote by 7c(M)] is the highest, but for 7 < the 
largest in module is the critical current 7c(^'^) of the vortex Consequently, in a field 
Hb = 1.7 the positive critical current of the junction is 7c = 7c(M), while the negative 
critical current is ■jc = id^^)- 
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Figure 5: The critical curves of the junction for the in-line (left) and overlap (right) 
geometries 

We note that for the in-line geometry the curve for the M distribution has a charac- 
teristic plateau — the "breakoff" of the Meissner solution sets in for a rather large module 
of the external current. 

The critical curve 7c(/i_B) for a junction is constructed as the envelope of the critical 
curves corresponding to different magnetic flux distributions in the junction. In other 
words, the critical curve consists of pieces of the critical curves for individual states with 
the largest module of the critical current for a given Hb- The parts of the critical curves 
corresponding to the intervals hs & [0, 4) for in-line and /i_b G [0, 5) for overlap geometries 
are illustrated in Fig. We note that our numerically constructed critical curve of 
the Josephson junction for the in-line geometry matches well with the theoretical and 
experimental results presented in Figs. 6a and 7a of Ref. S]. Critical curves in case of 
overlap ESJJ's have not been obtained experimentally yet. 

3 Conclusions 

Numerical calculations of the magnetic flux distributions and their bifurcation with pa- 
rameters of the model in the long ESJJ have been done. To study the stability of the 
solutions, we associate with each distribution the SLP whose potential is defined by this 
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distribution. We showed that the magnetic flux can change its stabihty for some critical 
values of the parameters. The different behaviour of the bifurcation upon a change in 
the shape parameter for different values of the boundary magnetic field and the external 
current was investigated. The envelope curve of the critical curves for several fiuxon states 
was obtained as the critical curve for the junction in cases of in-line and overlap geome- 
tries. Our numerically constructed critical curve of the exponentially shaped Josephson 
junction matches well with the experimental data. 
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